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Abstract. ATSM are widely applied for pricing of bonds and in- 
terest rate derivatives but the consistency of ATSM when the short 
rate, r, is unbounded from below remains essentially an open ques- 
tion. First, the standard approach to ATSM uses the Feynman-Kac 
theorem which is easily applicable only when r is bounded from 
below. Second, if the tuple of state variables belongs to the region 
where r is positive, the bond price should decrease in any state 
variable for which the corresponding coefficient in the formula for 
r is positive; the bond price should also decrease as the time to 
maturity increases. In the paper, sufficient conditions for the ap- 
plication of the Feynman-Kac formula, and monotonicity of the 
bond price are derived, for wide classes of affine term structure 
models in the pure diffusion case. Necessary conditions for the 
monotonicity are obtained as well. The results can be generalized 
for jump-diffusion processes. 
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1. Introduction 

Consider a financial market under several sources of uncertainty rep- 
resented by a multi-variate Markov process X. The price of an interest 
rate derivative of the European type, maturing at date T, with the 
terminal pay-off g(X(T)), can be expressed as 



f(X(t),t) = E t 



T 



exp - / r(X(s))ds g(X(T)) 



(1.1) 



Starting with one-factor diffusion models Vasicek (1977) and Cox, In- 
gersoll and Ross (1985), one of the popular approaches has been to 
model X as the solution to the stochastic differential equation 



dXj(t) = b J (X(t),t)dt + J2^S j (X(t),t)dW j} (1.2) 

k=l 

j = 1, . . . ,n, where bj and Sj are affine functions of X(t), and dW 
is the increment of the standard n-dimensional Brownian motion; r is 
modelled as an affine function of the state variable: 

r(X{s)) = (d,X(s)) + d , (1.3) 

where d G R™ is a constant vector, and d Q G R is a scalar; (•, •) 
denotes the standard inner product in R n . When r is given by ()1.3|) . the 
Feynman-Kac formula and the Fourier transform can be used to reduce 
the calculation of f(X(t),t) to the solution of a parabolic equation, and 
then to solution of a system of ODE (Riccati equations) depending on 
a parameter, with the subsequent integration w.r.t. this parameter. 
This idea is due to Heston (1993) who applied it to pricing of bond 
and currency options. Heston's approach was generalized by Duffie 
and Kan (1996), who coined the term Affine Term Structure Models 
(ATSM). For the classification of ATSM under diffusion processes, see 
Dai and Singleton (2000), and for the extension of ATSM to some jump- 
diffusion processes and extensive bibliography on different families of 
ATSM for both pure jump and jump- diffusion cases, see e.g. Duffie, 
Pan and Singleton (2000) and Chacko and Das (2002). Notice that the 
presence of jumps imposes additional restrictions on the parameters of 
the model. For instance, in one- dimensional case, one must ensure that 
jumps cannot move X(t) in the region where the volatility coefficient 
becomes negative. Thus, either the volatility is independent of the 
state variable, or an appropriate restriction on the direction of jumps 
must be imposed. For very general classes of affine Markov models 
with jumps, under conditions which ensure the non-negativity of r, see 
Duffie, Filipovic and Schachermayer (2002). 
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The consistency of ATSM in cases when r may be unbounded from 
below remains essentially an open question. The main stress in the 
classification paper Dai and Singleton (2000) is on the overdeterminacy 
of many ATSM models; however, for wide regions in the parameter's 
space, standard ATSM models may be inconsistent, and the following 
issues must be addressed. First, the standard approach is based on 
the Feynman-Kac formula but the general Feynman-Kac theorem is 
formulated for bounded (and sufficiently regular, say, continuous) r and 
sufficiently regular g; an extension to the case of r which are bounded 
from below is straightforward. Thus, the first step of the solution of 
an ATSM, namely, the reduction to the backward parabolic problem 

(d t + L-r)f(x,t) = 0, t<T; (1.4) 
f(x,T) = g(x), (1.5) 

where L is the infinitesimal generator of X, cannot be easily deduced 
from a general Feyman-Kac theorem unless the affine r depends on 
Xj of the CIR-type only 1 . If some of Xj may assume arbitrary (real) 
values, and the corresponding coefficients dj in (jl.3J) are non-zero, then 
r is unbounded from below, and to the best of our knowledge, no 
universal statement exists which guarantees that the solution to the 
problem (jl.4j) - ([1.5j) . call it fo(g, r; x, t), coincides with f(g, r; x, t) given 
by the stochastic expression (jl.lj) . Second, it is natural to assume 
that in the fully consistent model, the solution to the bond pricing 
problem must be a decreasing function of any state variable for which 
the corresponding coefficients in the formula for r are positive; the 
solution must also decrease as the time to maturity increases, if the 
tuple of state variables belongs to the region where r is positive. 

In the classification of Dai and Singleton (2000), a model is said to 
belong to family A m (n) if the number of the factors of the CIR-type is 
m. The two opposite cases, m = n and m = 0, are especially simple. If 
m = n, then the justification of the use of the Feynman-Kac theorem is 
a special case of the general result, and if in addition, d > 0, then the 
monotonicity of the bond price is evident from (jl.lj) . If in = 0, then X 
is the Ornstein-Uhlenbeck process on R n , for which an explicit formula 
for the characteristic function E x [e*^ ,x ^] is available (see equation 
(17.4) in Sato (1999) 2 ). The formula holds even for a wider class of pro- 
cesses of Ornstein-Uhlenbeck type, which are driven by Levy processes, 
and it is equivalent to the statement that the stochastic expression (jl.lj) 

1 We say that a factor Xj in an ATSM is of the CIR-type iff it assumes values in 
R+ 

2 The author is grateful to Darrell Dufhe for the reference 
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with g = 1 is equal to the solution to the Riccati equations for the bond 
price, that is, /o(l, r; x, t) = /(l, r; x, t). However, there is no reason to 
expect that the bond price is monotone in this case for all parameters' 
values. For 1 < m < n— 1, both the reduction to the Riccati equations 
and the monotonicity conditions have not been studied. 
Let r = T — t be the time to maturity, and 



be the price of the bond, which is obtained in an ATSM model by 
the formal reduction to the Riccati equations. First, we consider the 
Vasicek model and its generalization, namely, family A (n), then a 
simple two-factor Ai(2) model, next more general Ai(n) model, and 
finally the family ^(3) (other families A m (n) can be studied similarly), 
and derive, in terms of parameters of the model, 

(I) simple necessary conditions for the decay of P(x,t): 

P;(x,r)<0, Vx s.t. r(x)>0, and V r > 0; (1.7) 

in some cases, we also show that these conditions imply the 
boundedness of the bond price: in the region {x \ Xj > : 



(II) sufficient conditions for the decay of Bj(r); we do not know how 
wide is the gap between these conditions and the (unknown to 
us) necessary and sufficient conditions; 
(III) sufficient conditions under which the reduction to the system of 
the Riccati equations can be justified. For Ax(2) family, and in 
many other cases, these condition are weaker than the necessary 
conditions in (I). 

Remark 1.1. a) As it was mentioned above, for v4o(ri)-model, the re- 
duction to the system of Riccati equations is known, and it is valid 
without additional conditions on parameters of the model. 

b) Necessary and sufficient conditions for (jl.7|) in a neighborhood of 
r = +oo, and in a vicinity of 0, are easier to derive, and under these 
conditions, a "numerical proof of the monotonicity of the bond price 
on a large finite interval can be used to show that for given parameters' 
values, the model is consistent. 

c) As our study shows, for family Ai(n), the monotonicity of P(x, r) 
w.r.t. t is the main consistency problem for ATSM (and the only 
consistency problem for family Ao(n)). On the other hand, should one 
use the model for a fixed (and sufficiently small) time to maturity then 




(1.6) 



dj > 0}, 



P(x,t)<C, Vr>0; 



(1.8) 
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the model can be consistent on this time interval; and it is possible to 
derive sufficient conditions for (jl.7j) to be valid on a sufficiently narrow 
interval (0, r ), where r > depends on parameters of the model. 

d) When it is necessary to consider more general contingent claims, 
a sufficient condition for (III), in terms of the rate of growth of the 
pay-off at infinity, can be derived relatively easily, and the same is true 
of a necessary condition for the natural analog of (jl.7j) and sufficient 
condition for (jl.8|) . The sufficient conditions for the monotonicity will 
be more difficult to derive. 

e) It is plausible that in some empirical studies, the fitted ATSM is 
inconsistent in the sense that the monotonicity condition fails. Hence, 
if the model is fitted for some time to maturity, and used later for a 
larger time, then it may produce non-monotonic bond prices. 

f ) Similar consistency problems exist for interest rate derivative prod- 
ucts, and an unnatural behavior of the price of a derivative product can 
be easily overlooked if one fits the parameters of the model by using the 
data on bond prices, and then uses the calibrated model to calculate 
prices of interest rate derivatives. It might be possible to construct 
an arbitrage strategy against a counterparty who uses an inconsistent 
model. 

One may argue that the consistency analysis should be conducted 
when not only parameters of the model are fixed but values of (un- 
observable) factors as well: if for a chosen set of the parameters and 
factors the bond price is a decaying function of r, then the model is 
reasonable. However, it is not clear how to obtain a general result 
in this set-up, and moreover, it does not seem to be right to put the 
factors on the equal footing with the parameters of the model, espe- 
cially in cases when the factors can be interpreted as the short rate, its 
volatility and/or central tendency. 

We assume that the model should be consistent for all positive values 
of factors, the interest rate depends on (we consider the case when in 
(jl.3|) . the coefficients dj,j = 1, . . . , n, are non-negative), which makes 
the model less flexible than it is assumed in Dai and Singleton (2000). 
For the factors of the CIR-type, our restriction on g^ 's is natural, and it 
is without loss of generality as far as the factors assuming values in R 
are concerned: if one of these dj is negative, one can make it positive by 
using — Xj instead of Xj. On the other hand, the model specifications 
are underdetermined in the sense of Dai and Singleton (2000) because 
this allows us to formulate necessary conditions and sufficient ones in 
more symmetric and natural forms. 

The rest of the paper is organized as follows. In Section |2J we formu- 
late and prove the results for the Vasicek model and its generalization, 
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Ao(n)- model. Although the justification of the use of the Feyman-Kac 
theorem can be made by appealing to formula (17.4) in Sato (1999), 
we give an independent proof of the reduction. The main idea of the 
proof is the conjugation with an appropriate exponent, which allows 
us to reduce to the case of an interest rate bounded from below. The 
same idea is used in the proofs for A m (n)-models, m > 1, however 
the realization becomes much more involved, if we want to avoid too 
stringent conditions. 

In Section |2l we formulate the main results for the case m > 1, 
and prove theorems about properties of the formal solution (|1.6|) of 
the ATSM model. In Section we describe the general scheme of 
justification of the reduction to the backward parabolic equation; in 
other words, the scheme of the proof of the Feynman-Kac formula for 
an affine r and the class of processes used in A m (n) models, m > 1. The 
main ingredient of the proof is the representation theorem for analytic 
semigroups. The proof of the reduction for Ai(n) model and A 2 (3) 
model is given in Section (more general y4 m (ra)-models can be studied 
similarly). The most technical part of the proof, namely, the proof of 
the existence and uniqueness theorem for degenerate elliptic operators 
with parameter is delegated to the appendix. The proof of the latter 
theorem uses a general approach to problems of this sort described 
in detail in the review paper Levendorskii and Paneyakh (1990) and 
monograph Levendorskii (1993). This approach is applicable not to 
differential operators only but to integro-differential operators (another 
name: pseudo-differential operators) as well, which allows one to justify 
the use of the Feynman-Kac formula for jump-diffusion processes. This 
more general case will be treated in a separate publication. 



2. Family A (n) 

2.1. Monotonicity conditions. The interest rate is given by (jl.3|) 
with non- negative dj,j > 1, and positive d n , and the dynamics of X(t) 
is given by SDE 

dX(t) = (6-KX(t))dt + 5W(t), (2.1) 

where W(t) is the standard Brownian motion in R n , k = [kji] is a 
low-diagonal matrix with positive diagonal elements, S is a positive- 
definite matrix, and 9 is a constant vector. The bond price is given by 
(JEH with 

£?(r) = -(l-exp(-™ T )) 7 , (2.2) 



CONSISTENCY CONDITIONS FOR ATSM 



7 



where 7 = (k t ) 1 d, and 

C(r) = jT (~d + 6 T B(s) + ilr (E T [Bj(s)Bi(s)]^ l=1 T^j ds, (2.3) 

and necessary and sufficient conditions for the decay of the bond price 
are (relatively) easy to establish. These conditions are especially simple 
in the case n = 1 (the Vasicek model). 

Theorem 2.1. Let n = 1, and let the dynamics of the short rate be 
given by 

dr{t) = (8- nr{t))dt + adW(t), 

where 9, k and a are positive constants. 

Then a) if for some r > 0, the bond price P(r, r) is a non-increasing 
function of r, then 

y < (2-4) 

b) if \2.J$ holds, then for any r > 0, the bond price P(r, r) is a de- 
creasing function in t. 

Proof. Clearly, B(t) = — k~ 1 (1 — e~ KT ) is decreasing on [0, +00) from 
to — k~ 1 , therefore if 

C'(r) = 6B(t) + °^B{rf 

is non-positive on (0, +00), then (|2.4jl holds, and if ()2.4j) holds, then 
C'(t) is negative on (0, +00). □ 



Condition ()2.4j) means that for a given central tendency and coeffi- 
cient of mean reversion, the volatility may not be too large, and the 
interpretation is clear: if volatility is large, then a trajectory of the 
process spends significant amount of time in a region of the state space 
where the short rate is negative, which leads to the artificial increase 
of the bond price. 

In the case n > 2, it is convenient to study the monotonicity of 
Bj step by step, by using the low-diagonal structure of k, instead of 
appealing to the explicit general formula ()2.2|) . Due to positivity of the 
diagonal elements of matrix k and the assumption d n > 0, the B n is 
decreasing from to —d n /K nn . Hence, if 

Bn-i(r)= [ e-^-^-^-^i-K^B^s) - d n -x)ds 
Jo 

is non-increasing on (0, +00), then it must be that 

d n -l K nn ~ d n K n n _i > 0, (2-5) 
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and if (|2.5j) holds, then -B n _i decreases. For Bj,j < n — 2, the necessary- 
conditions for the monotonicity are not that simple but the induction 
and the same consideration as above show that if the off-diagonal en- 
tries of the matrix k are non-positive then all Bfs are decreasing func- 
tions. From (|2.3jl . it is evident that then the bond price is a decreasing 
function in r if and only if the function 

R+ 3 s h- F(s) := -d + (9,B(s)) + ^Tr (p r [B j (s)Bi(s)]j l=1 E) E R 

is non-positive. Clearly, necessary and sufficient conditions for non- 
positivity of F in terms of the parameters of the model cannot be 
simple, however, relatively simple necessary conditions and sufficient 
ones (the latter more stringent than the former) are easy to formulate. 
From (J2.2j) . we find that B(+oo) = —7 := — (K T )~ 1 d, and therefore for 
any j = 1, . . . , n, and s > 0, we have Bj(s) E (—Jj, 0). Introduce the 
quadratic polynomial 

Q(y) = -d + (9,y) + ^Tr (£ T [ % -^ i=1 £) . 

If Q(— 7) > 0, then for sufficiently large s, F(s) > 0, and therefore C 
grows in a neighborhood of +00; this gives a necessary condition 

Q(-7) < 0. (2.6) 

Since the quadratic form in the definition of Q(y) is positive-definite, 
Q is non-positive on U 1 := [—71, 0] x • ■ ■ x [— j n , 0] (which implies that 
the integrand in ()2.3|) is negative a.e., and C decreases) if and only if 
it is non-positive at each vertex of £/ 7 , that is, 

Q(y)<o Vj/e{-7i,o}x...x{- 7n ,o}. (2.7) 

We have obtained 

Theorem 2.2. a) Let n = 2, and let B\ he non-increasing. 
Then \2. 5)) holds. 

b) Let n>2, and let all Bj be non-increasing. Then \2. 6]) holds. 

c) Let n > 2, let the off-diagonal entries of k be non-positive, and 
let ||3) hold. 

Then the bond price is a decreasing function of (x, r) in the region 
x > 0, r > 0. 

2.2. Justification of the Feyman-Kac formula. We assume that 
S is non-degenerate. 

Theorem 2.3. Let g be a continuous function, which does not grow 
too rapidly at the infinity 

ln(l + \g(x)\) = o(||x|| 2 ), asx->+oo. (2.8) 
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Then the expressions M.l\) . hi. 2^1 and M.ty are the unique solution 
to the problem fll.^| )-( TT3j) in the class of continuous functions, which 
admit the bound uniformly in t G [0, T\. 

As it will be seen from the proof, condition on g can be relaxed, and 
the solution is unique in a wider class of functions. 

Proof. By making an affine transformation of the factors, we may as- 
sume that there exist Co > such that 

(kx, x) > cq\\x\ | 2 . (2.9) 

Denote by fo(g, r, ;x,t) a solution to the problem (jl.4|) - (jl.5j) . 

Lemma 2.4. a) Let r and g be continuous, let g satisfy and r 

satisfy 

r(x) = o(\ \x\ | 2 ), as x — > oo. (2.10) 

Then a solution to the problem fli.^| )-( fT~5)) in the class of continuous 
functions f(x,t), which admit the bound uniformly in t e [0,T], 
exists and it is unique. 

b) Let r%, r 2 , . . . , r^, ■ ■ ■, be a sequence of continuous functions which 
satisfies H2.1U\) uniformly in N, and converges pointwise to a function 
r . Then 

fo(g,r N ;x,t) -»• f (g,r;x,t) as iV -»• +00, (2.11) 

pointwise. 

Proof. The key element of our approach is the conjugation with an 
appropriate exponential function; in the case of A (n), this function is 
especially simple. Take a small e > 0, and set 

fe(g,r;x,t) = exp(-e||x|| 2 )/ (fl',r;x,t), 

g e (x) = exp(-e\\x\\ 2 )g(x). 

Similarly, define f e (g,r N ;x,t). Due to (J2HI), f e (g,r;x,t), f e (g,r N ;x,t) 
and g t {x) vanish as x — ► 00, uniformly in iV and t e [0, T]. 

a) Insert f (g,r; x,t) = exp(e| \x\ \ 2 )f £ (g, r; x, t) into (|1.4|) - (jl.5j) . and 
multiply by exp(— e| \x\ | 2 ). The result is a problem of the same form 

(d t + L e -r e )f £ (g,r,x,t) = 0, t < T; (2.12) 
fe(g,r,x,T) = g e (x), (2.13) 

where r e is a function, and L e is a differential operator without the first 
order term, which are obtained from 

exp(— e| \x\ \ 2 )(L — r) exp(e| \x\ | 2 ) = L t — r t . 

It is easily seen that the coefficients of L e at derivatives of order 2 are 
the same as the ones of L, and the other coefficients of L e tend to the 
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corresponding coefficients of L as e — > 0. Hence, if e > is sufficiently 
small, then L t is the infmitesial generator of an Ornstein-Uhlenbeck 
process, call it X e . Fix such an e. Further, r t = r + f e , where 

r e (x) = 2e{KX, x) + e(Tr(S T S) - 2(9, x)), 

and in view of (j2.1()j) and (|2.9jl . there exist positive c, C such that 

c||x|| 2 -C < r e (a;) < c| |a;| | 2 + C, VxGR". (2.14) 

Since r e is bounded from below, the solution to problem (|2.12|) - (|2.13l) 
exists in the class of continuous functions decaying as x — > oo, and it 
is unique (in fact, the solution is unique in the class of functions which 
grow not faster than an exponential function). Since fo(g,r; x,t) = 
exp(e| \x\ \ 2 )f t (g, t\ x, t), part a) has been proved. 

b) The argument above is applicable with instead of r, and 
rjv,e — r N + ?e admits bound (J2.14|) uniformly in N. It follows that 
fe{g, r N,x,t) is given by the Feynman-Kac formula 



fe(g,r N ,x,t) = E t 



exp(- / r N ^(s))ds)g e (X*(T)) 



(2.15) 



Since r^^x) — > r e (x), as N — > +oo, point- wise, we use the Domi- 
nant Convergence Theorem, pass to the limit in (|2.15|) . and obtain that 
f e (g, r N , x, t) -> f e (g, r, x, t), point-wise. It follows that f (g, r N , x, t) 
fo(g, r, x, t), point-wise. □ 

Now we give the proof of Theorem 12.31 Denote f(g,r;x,t) as the 
stochastic expression (jl.lj) . Fix N G R, and for the affine r, set 
tn{x) = max{r(x), — N}. Then rjy is bounded from below, hence both 
f(g,r N ;x,t) and f (g,r N ;x,t) exist, and f(g,r N ;x,t) = f (g,r N ;x,t) 
for all x and t < T. Since g, and r satisfy the conditions of 
Lemma 12.41 (j2.11j) holds. By the Monotone Convergence Theorem, 
f(g, r N ; x, t) f(g, r; x, t) as iV — >• +oo, point-wise, and we conclude 
that f(g, r; x, t) = fo(g, r; x, t) for all x and t <T. □ 

3. Families A m (n), 1 < m < n — 1. 

3.1. Family Ai(2). The state space is R + x R, the r is given by (jl.3|) 
with d\ > 0, c?2 > 0, and the infinitesimal generator of the process is of 
the form 

L = (6> x - KxiXx)di + (fl 2 ~ «2igi - K22X2)d 2 + + a + ^ Xl 8$, (3.1) 

where Kn, K22, 0i, ct, ft are positive. Without loss of generality, O2 = 0. 

Assume that in the case of the bond, the use of the Feynman-Kac 
theorem has been justified. Then P(x,t) = /o(l, r; x, t), the solution 
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to (jl.4|) - (jl.f)j) with g(x) = 1, can be found in the form (jl.6j) . By substi- 
tuting (jl.6|) into (jl.4)) - (jl.5j) . we obtain the system of Riccati equations 
on (0,T): 

B[ = - Klx B x + X -B\-K 21 B 2 + ^B\-d x , (3.2) 
B' 2 = -K 22 B 2 -d 2 , (3.3) 
C = -d + 9 1 B 1 + ^Bl (3.4) 

subject to the boundary conditions 

Bi(0) = 0; B 2 (0)=0; (7(0) = 0. (3.5) 
We solve Q subject to B 2 (0) = 0: 

J B 2 (r) = - 7 (l-e- K22T ), (3.6) 
where 7 := K 22 d 2 > 0. -B2 decreases from to —7, since K22 > 0: 

lim B 2 (t) = -7 < 0. (3.7) 

If Bi(tq) > for some tq > 0, then not only P(x,t) fails to be a 
decaying function of r; for this r , -P(x, r) is an increasing function in 
X\. The following theorem provides necessary conditions which exclude 
such a strange behavior of the bond price (the higher the spot short 
rate, the higher the price of the bond), and sufficient conditions for the 
negativity of B\(f) for r > 0. 

Theorem 3.1. a) If B\ is non-increasing on [0,+oo) ; then the follow- 
ing two conditions hold: 

d 1 > or k 2 i < 0; (3.8) 



di - k 2 i7 - 
b) If /TO)- /TO) hold, then 

B x {r) <0, 

Proof, a) Suppose that d\ = k 2 i — but Bi{r) is non-positive in a 
right neighborhood of 0. From ([3.6)1 . the RHS in ()3.2|) is positive in 
this neighborhood, hence B\ is increasing in this neighborhood from 0. 
Hence, B x is positive there; contradiction. Thus, (|3.8|) holds. 

Denote by d the LHS of flUSJ). In view of (jSZ3), the RHS of (Q 
admits the representation 

-KuBxir) + Ib^t) 2 -d + o(l), as r -> +00. 



P 2 

-7 2 > 0. 

2 ' ~ 



V r > 0. 



(3.9) 
(3.10) 



12 



S. LEVENDORSKrf 



Hence, if d is negative, and B\ stays negative, then eventually, the RHS 
of (|3.2|) will exceed a positive constant, and B\ will grow as a linear 
function of r, contradiction. 

b) From ()3.8|) and (|3.6|) . we see that Si is negative in a right neighbor- 
hood of 0. Consider the quadratic polynomial Q(y) = —di — K 2 iy + ^y 2 . 
Under condition (|3.8|) . it is negative in a small left neighborhood of 0, 
therefore if (|3.9|) holds, it is negative on (—7,0). Hence, 

-d x - k 21 B 2 (t) + ^B 2 (r) 2 < 0, V r > 0. 

We conclude that if -Bi(t) tries to approach zero from below as r in- 
creases, then the RHS in (|3.2j) becomes negative, and pushes the tra- 
jectory r 1 — > Bi(t) down. Hence, -Bi(t) remains negative. □ 

The sufficient conditions for the monotonicity are stronger than ()3.8j) - 
(13. 9|) ; in fact, additional necessary conditions can be derived. We are 
satisfied with ()3.8|) - ()3.9|) for the time being because these conditions 
are more restrictive than the sufficient conditions for the use of the 
Feynman-Kac formula in Theorem 13.31 below. The following theorem 
provides sufficient conditions for the decay of B\ and C . 

Theorem 3.2. a) Let \3. 8\) and \3. hold, and let 

/?7 + ^2i<0. (3.11) 

Then Bi decreases on [0,+oo). 

b) Let and KTTty hold, and let 

rv 

d > - 7 2 . (3.12) 

Then C decreases on [0, +00), and the bond price is a decreasing func- 
tion of (2, r) in the region Xi,x 2 > 0, r > 0. 

Proof, a) Set Y = B[, and solve ()3.2|) w.r.t. B 1 taking into account 
that the latter is negative: 

B x = «n - yjn 2 n + 2k 12 B 2 - (3B 2 + 2d x + 2Y. 

By differentiating (|3.2jl . we find 

Y' = -k 11 Y + YB 1 -k 21 B' 2 + PB 2 B 2 (3.13) 

= -Y^k 2 u + 2k 12 B 2 - PBI + 2c?i + 2Y + B' 2 ((3B 2 - « 2a ). 

On the strengh of ()3.8|) and (|3.6jh K is negative in a small right neigh- 
borhood of 0, and for all r > 0, B' 2 {r) < and B 2 {r) G (—7, 0). Hence, 
under condition (f3~TTj) . B' 2 (t)((3B 2 (t) - k 21 ) < 0, V r > 0, and if F(r) 
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tries to approach from below, Y' (r) becomes negative. Thus, B[ = Y 
must remain negative on (0, +oo). 

b) The RHS of (Q is less than -d + a 7 2 /2. □ 

The next theorem justifies the use of the formal solution of ATSM. 

Theorem 3.3. Let g be bounded and continuous, and let 

di + ^f- ^217 - f 7 2 > 0. (3.14) 

Then the expressions M.l\) . M.2]) and M.3]) are the unique solution 
to the problem fli.^| )-( f?~5j) in the class of continuous functions, which 
admit the bound 

|/(x,r)|<Cexp[(-7x 2 ) + ]. (3.15) 

Proof. The general scheme of the proof is described in Section HJ and 
the proof itself is given in Section |S| □ 

Remark 3.1. The bound (j3.15|) . which specifies the class of functions 
among which the Feynman-Kac formula gives the unique solution, is 
natural because the bond price given by (jl.6|) satisfies this bound: see 
(ESD. 

Remark 3.2. The LHS in the condition ()3.14|) summarizes the influence 
of different parameters of the process and the affine model of the short 
rate, which act in opposite directions. The main obstacle for the proof 
of the Feynman-Kac theorem (and the validity of this theorem) comes 
from the region in the state space where the short rate is unbounded 
from below. The larger the d\, the smaller this region is, and the 
condition (|3.14j) provides the lower bound for di, which ensures the 
validity of the Feynman-Kac (the other parameters remaining fixed). 
In the Aq(u) case, when the volatility is bounded and so trajectories 
of the process are not pushed too far too fast in the "negative" region, 
no additional condition is needed. In the Ai(n) case, the entries of the 
volatility matrix are large in the region where x\ is large, therefore the 
less time a trajectory of the process spends there, the smaller obstacle 
for the proof of the Feynman-Kac theorem (and the validity of this 
theorem) is. A large value of the mean-reverting entry Ku ensures that 
trajectories of the process do not remain far from the line X\ = for 
too long time, and the condition (|3.14j) provides the lower bound for 
Ku, the other parameters remaining fixed. 

The parameters <i 2 , P and k 22 act in opposite directions as well. The 
larger the cf 2 , the larger (the absolute value of) negative values of the 
short rate can be, and the larger the f3, the farther a trajectory can be 
pushed into the "negative" region. The larger the k 22 is, the faster a 
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trajectory of the process will return to a neighborhood of a line x 2 = 0, 
where the short rate is bounded from below. If k 2 i — 0, then (J3.14j) 
provides an upper bound on [3{d 2 / K22) 2 , which ensures the validity of 
the Feynman-Kac theorem. If k 2 i 7^ 0, the dynamics becomes more 
complex, and the term K217 takes this additional complexity into ac- 
count. 

Corollary 3.4. For the bond pricing problem, Theorem VSJA is valid 
under condition 

d i + ^f- K 2i7 - ^7 2 > 0. (3.16) 

Proof. Indeed, suppose that (j3.16J) holds with the equality. Then we 
can approximate di by a sequence {d lm } converging to di from above. 
Each di <m satisfies (J3.14j) . and hence the conclusion of Theorem 13.31 
holds. Set r m (x) = d + d lm x\ + d 2 x 2 . By the Dominant Convergence 
Theorem, f(r m ,g;x,t) — > f(r,g;x,t), point-wise, and the convergence 
fo( r m, g', x, t) — >■ f Q (r, g; x, t) follows from the explicit formula for these 
functions and the theorem about the continuous dependence on param- 
eters of the solution of a system of ODE of the first order. □ 

3.2. Family Ai(n). The state space is R + x R n_1 , the r is given by 
(jl.3j) with dj > 0,j = l,...,n, and d n > 0, and the infinitesimal 
generator of the process is of the form 

1 1 n 

L=(B- nx, d) + -x x dl + 2 ( a i k + x iPjk)djd k , (3.17) 

j,k=2 

where 9\,Kjj,j = 1, . . . ,n are positive, Kji < 0, 1 < I < j, Kjk — 0, j < 
k, and a = [otjk]"j t k=2-> & = [Pjk) 7 j,k=2 are positive definite matrices (these 
restrictions can be relaxed). Without loss of generality, we may assume 
6j = 0,j > 2, when convenient. 

Introduce vectors d 2 = (d 2 , . . . , d n ), k 21 = (k 2 i, ■ ■ ■ , ^ni), and set 7 = 
(K T )~ 1 d. The subsystem of the Riccati equations for B = (B 2 , . . . , B n ) 
enjoys the same properties as the subsystem for B in the A (n)-model 
above, therefore the same proof as of Theorem 12.21 gives 

Lemma 3.5. For j = 2,...,n, function Bj is continuous and non- 
increasing on [0, +oo), and Bj(0) = 0,Bj(+oo) = —jj. 

It remains to study the monotonicity of B\ and C . Introduce qua- 
dratic polynomials 

Qb(v) = ~dt - (y, k 21 ) + ±y T [3y, y e R n ~\ 
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and 



Qc{y) = -do + ^y T ay, y G R n l . 



The remaining two Riccati equations, for B\ and C, are 

1 



k 11 B 1 (t) + ~B 2 (t) + Q b (B(t)), 



(3.18) 
(3.19) 



B'(r) = 

C\t) = 9 1 B 1 (t) + Q c (B(t)). 

Each of equations ()3.18j) and ()3.19j) is similar to equation ()2.3j) for C 
in the Ao(n)-model, therefore the same proof as of Theorem 12.21 gives 

Theorem 3.6. a) If B\ is non-increasing on [0,+oo) ; then 

di>0 or (d 2 ,K 21 ) < 0, (3.20) 

and 

Q B (- 72 ,...,- T „) < 0. (3.21) 

b) Let \3.2U\) hold, and 

Q B {y) < 0, V y E {-72, 0} x ■ • • x {- 7n , 0}. (3.22) 

Then B\ decreases on [0, +oo). 
c) Let pPfy) and 0JZ$) hold, and 

Qc{y) < 0, V y E {-72, 0} x ■ ■ • x {- 7n , 0}. (3.23) 

Then C decreases on [0, +oo) ; and the bond price is a decreasing func- 
tion of (x, r) in the region x > 0, r > 0. 

By using the family of transformations described in Dai and Single- 
ton (2000), it is possible to reduce any Ai (n) -mo del to a model with n 
satisfying 



K jk = 0, Vk^j, k>2, j> 1, 



(3.24) 



and we will formulate the theorem about the justification of the use of 
the Feynman-Kac theorem under this assumption. 

Theorem 3.7. Let g be bounded and continuous, let \3.2J^ ) hold, and 
let 

^f-Q B (y)>0, V y G {-72, 0} x • • • x {- 7n , 0}. (3.25) 

Then the expressions M.l\) . M.2]) and M.3]) are the unique solution 
to the problem Ji.^| )-( f?~5)) i n the class of continuous functions, which 
admit the bound 



\f(x,r)\<Cexp 
Proof. In Subsection 15.31 



J=2 



(3.26) 
□ 
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Remark 3.3. The bound (|3.2fij) . which specifies the class of functions 
among which the Feynman-Kac formula gives the unique solution, is 
natural because the bond price given by (J1.6J) satisfies this bound: see 
Lemma 13.51 

Remark 3.4- In order to apply Theorem 13.71 to a particular model, the 
family of transformations described in Dai and Singleton (2000) should 
be used to obtain a model which satisfies (|3.24|) . The condition ()3.25|) 
is a natural generalization of the condition (j3.14j) for the multi-factor 
case; the interpretation is essentially the same as in Remark 3.2. 

Corollary 3.8. For the bond pricing problem, condition \3. 25}) in The- 
orem \3. 7\ can be replaced by a weaker one: 

^f-QB(y)>0, \/ye {-72, 0} x • ■ ■ x {- 7n , 0}. (3.27) 

3.3. Family A 2 (3). By using the family of transformations described 
in Dai and Singleton (2000), it is possible to reduce any three-factor 
ATSM with two factors of the CIR-type to a model of the form 



dX(t) = (6- KX{t))dt + W[Sy(t)]dB(t), (3.28) 

where 9 G R 3 is the vector with components #3 = 0, 

6 1 > 0, 9 2 > 0, (3.29) 
and the entries of matrix k = [kji] satisfy 

(3.30) 
(3.31) 
(3.32) 
(3.33) 



Further, 



where 



Kll, K 22 , ^33 


> 


0. 




< 


0. 


K 11 K 22 — K 12^21 


> 


0. 


= ^23 




0. 




J = 


1,2 



S 3J (t) = faXjit), J = 1,2, (3.34) 
S 33 (t) = a 3 + p 31 X 1 (t) + p 32 X 2 (t), (3.35) 



a 3 ,Ai,/?22 > 0, /3 3 i,/3 32 > 0, (3.36) 
and finally, the short rate process is given by (jl.3|) with 

d G R, di,d 2 >0, d 3 >0. (3.37) 

Notice that under conditions ()3.29|) and (j3.30j) - ()3.33p . any trajectory 
of the process X, which starts in the region X\ > 0, x 2 > 0, remains 
in this region, a.s. 
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The infinitesimal generator of the process is 

L=(6- Kx) T d x + - ^ hj x j d *j + ^("3 + An^i + (3 32 x 2 )dl, 

i=i,2 

and therefore the system of Riccati equations is 

B[ = -k 11 B 1 - k 2X B 2 + f^Bf -dx- k 31 B 3 + SEBl (3.38) 

B' 2 = -k 12 B 1 - k 22 B 2 + ^-B\ -d 2 - k 32 B 3 + ^Bl (3.39) 



B' 3 = (3.40) 

C = -d + 6 1 B 1 + 6 2 B 2 + ^B 2 3 . (3.41) 

From the initial condition 

^(0) = B 2 (0) = B 3 (0) = C(0) = (3.42) 
and (J3.4U)) we can easily find B 3 ; 

B 3 (r) = - 1 (l-e-™ T )i (3-43) 
where 7 = d 3 / k 33 > 0. Clearly, B 3 decays on [0, +00), and 

5 3 (+oo) = -7, (3.44) 

B 3 (r) e (-7,0), Vr>0, (3.45) 

B 3 (t) ~ -d 3 r, as t -* +0. (3.46) 



The monotonicity of Bj, j = 1, 2, implies additional restrictions on the 
parameters of the model. To formulate them, denote by k 11 the inverse 
to the upper left 2x2 block k 11 := [Kji]j i=i 2 of matrix k, 

^22 ~ K 12 
-K 2 i K n 



and set 



di 




d\ 




«31 




«31 




" hi ' 


= k 11 




. d-2 _ 




. ^2 . 


5 


. ^ 32 . 




. K32 . 


5 


. ^2 _ 




J32 _ 



On the strength of ()3.30|) - (j3.32|) . the diagonal (resp., off-diagonal) en- 
tries of k 11 are positive (resp., non-positive), therefore 

hi > 0, P32 > 0, (3.47) 



and 

di,d 2 >0 di,d 2 >0, (3.48) 
k 3 i,K32<0 =$> k 31 ,k 32 <0. (3.49) 



k n 



K 11 K 22 — K 12 K 21 
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Theorem 3.9. Let Bj,j = 1,2, be non-increasing on [0, +oo) 
the following conditions hold, for j ^ I G {1,2}: 

dj > or dj = and fyjdi + K3jd 3 < 0; 
dj > or = and k^j < 0; 



dj 



«3i7 - -y 7 



> 



0. 



Then 

(3.50) 
(3.51) 

(3.52) 



Proo/. a) If dj < 0, then (l3~3KJ) - (l3~39l and 1EQ5I) imply that Sj(r) > 
in a right neighborhood of 0, contradiction. Hence, both dj > 0. If 
both d j = 0, then from (13~3%1) - (13~3H|) and IfOBjl . both 5 3 -(r) = o(r), as 
r — ► +0, and moreover, 



-0. 



Hence, it is necessary that n^j < 0. Finally, if d\ = but g?2 > 0, then 
-82(7") ~ —d 2 r, as r — > +0, hence, 

fii(r) ~ («2id 2 + K3iG?3)t, r -> +0. 

This excludes the case ^21^2 + ^31^3 > 0, and finishes the proof of 

(EH) . 

b) To prove (I33TJ) - (J332|) . we apply /t 11 to subsystem (13~3%J) - (|3~331) : 

the result is 



B' 



B x 
B 2 




d x ' 
d 2 


-B 3 


«31 
^32 



(3.53) 



-Bi 



foi 
P32 



where 



B 



fri 
P22 



Since B[ and B' 2 are non-positive, and the entries of k u are non- 
negative, the LHS in (j3.53j) is non-positive. Now by using ()3.53|) and 
arguing as in part a) above, we deduce (|3.51|) . 

If (|3.52|) fails, then from (|3.44j) we conclude that in a neighborhood 
of +00, one of the components of the RHS in ()3.53j) is positive, contra- 
diction. □ 



The next theorem gives sufficient conditions for the monotonicity of 
Bj and C; they are more stringent than the necessary conditions in 
Theorem 13.91 
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Theorem 3.10. a) Let the following conditions hold, for j = 1,2: 

dj > or dj = and K 3 j < 0; (3.54) 

dj - K 3jl - ^ 7 2 > 0; (3.55) 

^ + Psj-y < 0, (3.56) 
fa > or K 3j < 0. (3.57) 

Then B\ and B 2 are decreasing on [0, +00). 
b) In addition, let 

do > y7 2 - (3.58) 

Then C decreases on [0, +00) ; and the bond price is a decreasing func- 
tion of (x, r) in the region x\ > 0, x 2 > 0, r > 0. 

Proof, a) First, we show that there exists r > such that for j = 1, 2, 

Bj(T) < 0, < r < t , (3.59) 

and 

Bj(r) < 0, < r < r . (3.60) 

We use (gIH - (El , (B32D and fl3~gjp . If d 1 and d 2 are positive, then 
B'-{t) ~ — dj, as r — > +0, hence (|3.60p hold, and ()3.59|) holds as well. 
If both dj = but both n 3 j < 0, then Bj(t) ~ n 3 jd 3 r, as r — > +0, and 
(13.59}) and ()3.6()|1 hold. Finally, if one of dj, say, d 2 , is positive, and the 
other, di, is 0, then B' 2 {t) ~ — d 2 , B 2 {r) ~ —d 2 r, asr-> +0, and 

Bi(r) ~ (k 2 i^2 + «3i4)t, t- -> +0. 

Since k 2 i < 0, g? 2 > and d 3 > 0, we obtain ()3.59|) - ()3.60|) . 

Second, we show that (|3.59|) holds with tq = +00. Under condition 
f!3.54j) . the polynomials 

Qj{y) = -dj - n 3j y+^-y 2 

are negative on (—7, 0), therefore in view of ()3.45|) . 

-d, + k 3i B 3 (t) + ^B 3 (t) 2 <0, Vr>0. (3.61) 

Suppose that as r increases, Bj(r),j = 1,2, start to approach from 
below, simultaneously. Then from (J3.38j) - (|3.39j) and (J3.61|) . at least one 
of B'jij) becomes negative before both Bj(r),j = 1, 2, reach 0, contra- 
diction. If Bi{r) is approaching but B 2 {r) is not, then eventually, on 
the strength of (|3.38|) and condition k 21 < 0, B[(r) becomes negative, 
contradiction. Thus, ()3.59|) holds on the whole half-axis. 
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It remains to prove that Yj := Bj are negative on the whole half- axis, 
j = 1,2. We differentiate 

Y( = -K 11 Y 1 -K 21 Y 2 + i3 11 B 1 Y 1 -K 31 Y 3 + i3 31 B 3 Y 3 , (3.62) 
F 2 ' = -k 13 Y 1 - K22 Y 2 + P 22 B 2 Y 2 -k 32 Y 3 + P 32 B 3 Y 3 , (3.63) 

and rewrite ()3.62j) - (j3.63|) as 

Y{ = -{k xx -(3 xx B x )Y x -k 2X Y 2 + {(3 3X B 3 -k 3X )Y 3 , 
Yi = -k X2 Y x -(k 22 -(3 22 B 2 )Y 2 + ((3 32 B 3 -k 32 )Y 3 . 

Since Bj are negative, Kjj — PjjBj > 0,j = 1,2, and due to (|3.45|) . 
()3.57j) and ()3.56j) . all the expressions in the brackets are positive. Since 
Y 3 is negative, Y{(t) and Y^ij) cannot approach from below simulta- 
neously: indeed, then eventually, the RHS's will become negative, con- 
tradiction. Since — k 2X > 0, we have — k 2X Y 2 (t) < where Y 2 (t) < 0, 
therefore by the same reasoning, Y x (t) cannot approach from below 
while Y 2 (t) remains separated from 0. By interchanging the indices 1 
and 2, we conclude that (J3.60|) holds on the whole half- axis. 

b) Recall that 6 X and 9 2 are positive, and apply (|3.41|) . (|3.45|) and 

(£23). □ 

Theorem 3.11. Let g be bounded and continuous, and let the following 
conditions hold: for y = 0, 7, 

d. + ^ + ^-^y-^y 2 > 0, (3.64) 

k\ 2 k X2 k xx (3 32 2 

d 2 + —- + — n 32 y-—y > 0. (3.65) 

Then the expressions M.l\) . Al.fy) and M.'J\i are the unique solution 
to the problem jl.^| )- P~3j) in the class of continuous functions, which 
admit the bound 

\f{x,r)\<Cexp[(-jx 3 ) + }. (3.66) 
Proof. In Subsection 15.41 □ 

Remark 3.5. The bound (j3.66|) . which specifies the class of functions 
among which the Feynman-Kac formula gives the unique solution, is 
natural because the bond price given by (jl.6J) satisfies this bound: see 
(|3~23l . 

Remark 3.6. The pair of conditions ()3.64|) - (j3.65|) is a natural general- 
ization of the condition ()3.14|) in the A x {2) model; the interpretation 
is essentially the same as in Remark 3.2. 
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4. Justification of the use of the Feynman-Kac formula 

4.1. General scheme. Step 1. We assume that g is non-negative, 
bounded and sufficiently regular, so that in the case of a continu- 
ous r bounded from below (not necessarily affine), problem (jl.4j) - (jl.5|) 
has a solution, f (g,r; ■, ■), in the class of bounded continuous func- 
tions, which is given by the stochastic expression (jl.lj) : fo(g,r; x,t) = 
f(g, r; x, t) for all x and t < T. For affine diffusions, the coefficients of 
the infinitesimal generator of the process satisfy the global Lipschitz 
condition, and hence, the Feyman-Kac theorem is applicable. (The 
reduction to the case of more general g is fairly standard). 

Step 2. Fix iV G R, and for the affine r, set r^{x) = max{r(x), —N}. 
Then is bounded from below, hence both f(g, r^; x, t) and fo(g, r^; x, t) 
exist, and f(g, r N ; x, t) = fo(g, r N ; x, t) for all x and t < T. 

Step 3. By the Monotone Convergence Theorem, 

f(g,r N ;x,t)->f(g,r;x,t) as iV -> +oo, 
point-wise, therefore it remains to show that 

f (g, r N ; x, t) -> f (g, r; x, t) as iV -> +oo, (4.1) 

point-wise. 

Step 4- We prove ()4.1|) by reducing to the case of a family of short 
rates, which is bounded from below uniformly in N. To this end, we 
take a non- negative function G C°°(R n ), and consider the represen- 
tation 

exp(-0(x))Lexp(0(x)) = L <i> - r^x), (4.2) 

where is a differential operator without the zero-order term, and 
is a function. In the presence of factors of the CIR-type, the infini- 
tesimal generator may have non-trivial affine coefficients at the second 
order derivatives, therefore the choice of </> in the form of a quadratic 
polynomial, as in the proof of Theorem 12.31 is impossible except for 
some very special cases. One may try 0's which are (approximately) 
positive homogenous of degree 1 in a neighborhood of infinity. The 
simplest version <p(x) = (Ax, x) / (x) , where (x) := (1 + ||x|| 2 ) 1//2 , is 
possible but it requires unnecessary strong restrictions on parameters 
of the model. It turns out that the construction of <p should be adjusted 
to each model. For the proof to work, the following general properties 
of cf) are essential (they can be relaxed, though): 

(i) there exist constants Co > and C such that 

<f>{x) > cq\x\ - C; (4.3) 
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(ii) there exist constants Co > 0, M, C\ such that := r + 

admits the following estimate: 

C \x\ - M < Too^x) < d\x\ + M, Vi; (4.4) 

(hi) there exist constants C2 and 5 > such that for all x and t G 
[0, T], the function f^(g,r;x,t) := exp(—<fi(x))fo(g,r;x,t) satis- 
fies an estimate 

\U(g,r,x,t)\<C 2 exp(-6\x\); (4.5) 

(iv) is sufficiently regular in the sense that for any r° which admits 
the bound (|4.4jl . and a continuous <?°(x) admitting the bound 

\g°{x)\<C 3 exp(-5\x\), (4.6) 

where 5 > and C 3 are independent of x, a continuous solution 
to the problem 

(d t + L^-r°)f(x,t) = 0, 0<t<T; (4.7) 
f(x,T) = g°(x), (4.8) 

which exponentially decays at infinity, exists and it is unique; call 
it f^ig , r°; x, t) (we add indices 0, and use labels r° and g°) in 
order to avoid the confusion with g and r in (|1.4|) - (|1.5Jl ): 

(v) as iV — > +00, 

fig , r N ^ x, t) -> f+(g°, r^- x, t), (4.9) 

point- wise; here '■= + ?>• 

Notice that the construction of will depend on a small parameter 
e G (0, 1), and the will satisfy a weak version of the global Lipschitz 
condition, with parameter, which simplifies the proof of the existence 
and uniqueness theorem in Section 03 We do not specify this condition 
here because it is a useful technical tool only, and it should be possible 
to prove the existence and uniqueness theorem under weaker conditions. 

When the existence and uniqueness theorem is proved, the conver- 
gence (J4.9)) is typically not difficult to establish. For instance, if it can 
be shown that L$ is the generator of a Markov process without killing, 
and the Feynman-Kac theorem is applicable, then (J4.9|) can be easily 
deduced from the Feynman-Kac formula and the Dominant Conver- 
gence Theorem. However, for the L^, which we will construct below, 
there is no ready Feynman-Kac theorem available, and so ()4.9|) will 
be proved differently, by using the representation theorem for analytic 
semigroups. 
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When the is constructed, and properties (i)-(v) are established, 
we can prove (|4.1|) as follows. For N = 1, 2, . . . , oo, set 

g^x) = exp(-</>(x))g(x), 
U(g,r N ;x,t) = exp(-(f)(x))f (g } r N ;x } t), 

substitute f (g,r N ;x,t) = exp((f>(x)) f^g , r N ; x, t) into ()1.4|) - (jl.5j) . and 
multiply by exp(— <p(x)); we obtain that f,p(g, Tn', ■, •) is a continuous 
solution to the problem 

(d t + L^-r Ni(P )f(x,t) = 0, 0<t<T; (4.10) 
f(x,T) = g^x), (4.11) 

which exponentially decays at infinity. Clearly, r^,^ satisfy ()4.4|) with 
the constants independent of N — 1, 2, . . . , oo; the RHS in (j4.8J) admits 
bound ()4.6|) since g is bounded, and <p satisfies ()4.3|) . By (iv), the con- 
tinuous solution to problem (|4.1Uj) - (j4.11|) . which exponentially decays 
at infinity, is unique, hence U{g,r N ; -, •) = r^; -, •), and by (v), 

as N -> +oo, /*(s>, r^; •, •) converges to /*(^, r^; •, •), pointwise; 
(jOl follows. 

4.2. Outline of the verification of conditions (iv)— (v). Let x' be 

the group of the CIR-type variables, and x" the group of the other 
variables; then x = (x',x"). Let (R + ) m x R™-" 1 be the corresponding 
decomposition of the state space. We apply the standard approach: 
first, we write problem (J4.1(J|) - (J4.11)) in the form 

F'(t) + AF(t) = 0, r>0; (4.12) 
F(+0) = g°, (4.13) 

where A = —L^ + r° is a (partial) differential operator on (R+) m X 
R n_m , and for each r, F(r) is a function on (R + ) m x R n ~ m . Next, by 
using the Laplace transform w.r.t. r, we reduce problem (|4.1U|) - (|4.11j) 
to the family of problems 

(A + A)F(X) = g°, (4.14) 

where F is the Laplace transform of F, and A belongs to a half-plane 
of the complex plane, of the form {A | 3A > Ao}. Then, by using the 
theory of degenerate elliptic operators with parameter, we show that if 
Ao is sufficiently large, then the operator A + A is invertible uniformly 
w.r.t. A in the half-plane 3A > Ao- This proves the existence and 
uniqueness of the solution to the boundary problem (|4.10|) - ([4.11|) (that 
is, part (iv) of the general scheme) but in a wider class of functions 
since under this approach, ()4.13|) is satisfied in a weak sense (in the 
sense of the theory of generalized functions). To show that the solution 
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satisfies (I4.1HJ) in the strong sense, and to prove the convergence in part 
(v), some additional effort is needed. We show that the representation 
theorem for analytic semigroups can be applied to A, and derive (|4.9|) 
from this representation. We need the following definition and theorem. 

For a G (0,tt) and A > 0, set Scr,Ao = {A G C | |A| > A , argA G 
[— a, a}}. Let A be an operator in the Banach space B, and let there 
exist cr g (0, 7r), Ao > 0, and C\ such that for A G £<t,a j A + A is 
invertible, and the resolvent satisfies the estimate 

||(A + A)" 1 !! < d(l + IAD" 1 . (4.15) 

Then A is called a weakly a-positive operator 3 . If C = 0, A is called 
a-positive. Let C a> \ = <9X CT Xo be a regular contour, with a parameter- 
ization A = X(t) satisfying argA(t) = ±cr for t in a neighborhood of 
±oo. 

The following theorem is a special case of the representation theorem 
for the analytic semigroups (see Section IX. 10 in Yosida (1964)). 

Theorem 4.1. Let A be an unbounded operator in the Banach space B, 
and let there exist a G (tt/2,tt) and C such that A is weakly a-positive. 
Then 

a) for any g° G B, and any t > 0, the following integral is well-defined 

exp(-rA)g = (2m)- 1 f e rX {\ + A)' 1 g Q d\- (4.16) 

b) \^ . 1 6| ) defines a strongly continuous semigroup {T T } T >o in BbyTo = 
I, T T = exp(— rA), t > 0, and 

c) for any g° G B, F(r) = exp(— rA)g° is a strongly continuous solu- 
tion to the problem 

As B, we take L2((R+) m x R n_m ), and we show that the conditions 
of Theorem 14. II are satisfied. For a function h, denote by h(-) or simply 
h the-multiplication-by-/i-operator. We will be able to prove that if 
p > is sufficiently small then 

||exp(p(r)(A + r (-))(A + A)- 1 exp(-p(.)")||<^, (4.17) 

for all A G S CT) a . By using the representation ()4.16|) and estimate 
(|4.17|) . we can prove (|4.9|) as follows. Set r° = r^^, A = —L^ + r^^, 
and let F be the solution to (|4.12|) - (j4.13j) . Let Fi N ^ be the solution of 
the problem (j4~T^) - (j4~T^|) with A {N) = —L^ + r N ^ instead of A. We 
show that A and each A^ satisfies (J4.15|) uniformly in N, with the 



Usually, the label 9 is used instead of cr; unfortunately, 8 is already occupied as 
the standard notation in ATSM's 
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same £<t,Ao, and therefore the representation (|4.1fi|) holds for F^) with 
A(at) instead of A. Notice that 

(A + A (JV) )-V-(A + ArV = (A + A w )- 1 (r OOi0 -r iV , )(A + A)-V 

= (\ + Ayr i )- 1 b N)X C x g 1 , 

where 

&at,a := _ r iV,4>) ex P(-p(-)")(-^ + r oo,</>)~\ 

C7 A = exp(p(-)")(A + W-))(A + A)" 1 exp(-p(->"), 
and g 1 := exp(p(-)")g° is a continuous function which decays at infinity, 
if p < 5. It is easily seen that 

lim sup \bjy \(x)\ = 0, 

therefore the norm of the-multiplication-by-fr^A operator tends to zero 
as N — > +oo, uniformly in A G £ a ,\ i the norm of C\ is uniformly 
bounded w.r.t. A G £<r,A 011 the strength of (|4.17|) . Hence, 

F (n) (t)-F(t) = exp(-TA (N) )g° - exp(-rA)g° 

= (2m)- 1 I e^ A ((A + A^)" 1 - (A + Ay 1 ) g°d\ 

= (2m)- 1 I e^X + A^y'b^xCxg'dX 

vanishes as iV — > +oo, and ()4.9|) is proved. 

The regularity theorem for locally elliptic operators guarantees that 
the solutions f(N)(x,r) = F^)(x) are not only in L 2 ((R+) m x R"~ m ) 
but continuous as well, and they decay at the infinity. Hence, the con- 
vergence in C(R + ; L 2 ((R+) m x R"~ m ), the space of continuous vector- 
functions with values in L2((R+) m x R n_m ), implies the pointwise con- 
vergence of f( N )(x, t) to f(oo)(x, t). 

In Sectional the construction of (ft and the proof of the existence and 
uniqueness theorem and bound (j4.17|) will be provided for the family 
Ai(2), and then for families Ai(n), n > 2, when the construction of 
(f) becomes more involved. In Subsection 15.41 the modifications of the 
proof for the family A 2 (3) are outlined, and the proof for other families 
A n (m) is essentially the same. 

5. Proofs of sufficient conditions for the Feynman-Kac 

THEOREM 

5.1. Proof of Theorem 13.31 part I: construction of cf> and veri- 
fication of conditions (J4.2|) — (14.51) . Fix any positive 7+ and 

7- < -7- (5.1) 
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Next, take a non- decreasing function x £ C°°(R) such that 

x{y) = i-i v<~ l \ 



(5.2) 



x(y) = i+, y> 0; 



(5.3) 



and 



7- < x(y) < 7+> Vi/eR. 



(5.4) 



Then, for any e G (0, 1), construct functions Xe(y) — x{ e u) an d 



and finally, set 

<f) e (x) = K n xi + ifi £ (x 2 ). 
Clearly, (fi e satisfies (14. 3 j) . Set 5 := max{7 + , Kn, — 7 — 7_}; it is positive 
since k%\ > 0,7+ > and (|5.1j) holds. 

Lemma 5.1. For any e G (0, 1), tnere exists C e such that for all X\ > 
and £2 G R ; 



Proof. From the explicit solution to the bond pricing problem, we 
know that f (r, l;x, t) admits the bound ()3.15|) . Since g is bounded, 
f (r,g;x,t) admits the same bounds (with a different C, perhaps). 
Hence, on the set {x \ X\ > 0,x 2 > — 1/e}, function f (r,g;-,-) is 
bounded, and since <j) e is bounded from below by a linear function with 
positive coefficients, KnXi + 7+^2, we conclude that estimate ()5.5j) 
holds on this set. On the set {x \ x\ > 0,x 2 < —1/e}, function 
f (r, g\ -, ■) is bounded from above by an exponential function of the 
form C exp(— 7x2), and <j) £ is bounded from below by an affine function 
^11^1 + 7-^2 + c, where 



In view of our choice (|5.1jl . estimate (|5.5jl holds on this set as well. □ 

Now we can check that for any sufficiently small positive 7+ and 
e G (0, 1), and 7_ in a sufficiently small left vicinity of —7, the function 
(fi t satisfies conditions ()4.2|) and ()4.4|) . The dependence on e will also be 
used to check conditions (iv)-(v) of the general scheme: it is convenient 
to use the dependence on a small parameter, and prove the existence 
and uniqueness of the solution of the boundary problem ()4.7|) - (j4.8J) for 
sufficiently small e > 0. 




exp(-<j> e (x 2 ))\f (r,g;x,t)\ < C e exp(-5(x 1 + |jc 2 |)). 
In particular, \4-£\) holds. 



(5.5) 
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We have 

exp(-0 e (x))Lexp(0 e (a;)) = (9 1 - K 11 x 1 ){d 1 + «u) 

-{k 2 \Xi + K 22 X 2 )(d 2 + Xe{x 2 )) 

+\x l (d 1 + Kll ) 2 + ^^(d 2 + X e(x 2 )f 

therefore 

= e^ + ^xdl (5.6) 



and 



2 

+ [-k 21 x 1 - k 22 x 2 + (a + (3xi)xe(x 2 )]d 2 + a + ^ Xl d 2 2 , 
a 



TooA X ) = d Q - OxKxi - ^(ex' e {.X 2 ) + Xe{.X2f) (5.7) 

k 2 8 

+ [dl + ~Y + K2lXe(x 2 ) - -(ex' e (x 2 ) + Xe{x 2 f )]Xi 
+ [d 2 + K 22 Xe(x 2 )]x 2 . 

Here x'Av) ■= (x')e(y) ■= 

Lemma 5.2. Let \3. 14\ ) hold. Then there exists eo > 0, 7+ > and 
7_ < —7 such that for alle G (0, €q), estimate {j4-4\ ) holds with constants 
Co,Ci independent of e G (0,eo), and M = M e _1 ; where M is also 
independent of e. 

Proof. In view of ()5.4j) . the first three terms on the RHS of ()5.7j) are 
bounded, and 

sup |e^(z 2 )| -> 0, e->0, 

therefore it suffices to show that 7+ > and 7_ < —7 can be chosen 
so that 

1) d\ + k\xI '2 + n 2 iXe(x 2 ) — Pxe(x 2 ) 2 /2 is positive and bounded away 
from zero, uniformly in x 2 and e > 0; 

2) d 2 + K 22 Xe{x 2 ) is positive and bounded away from 0, uniformly in 
x 2 > and e > 0; 

3) d 2 + n 22 Xt{x 2 ) is negative and bounded away from 0, uniformly in 
x 2 < — 1/e and e > 0. 

Now, due to (|5.2|) . 3) is d 2 + k 2 27- < 0, which is equivalent to (|5.1|) . 
and due to (|5.3|) . 2) reduces to <i 2 + ^227+ > 0, which is trivial. On the 
strength of ()5.4|) . 1) is equivalent to the condition: the polynomial 

V >-> <5(y) := di + ^ + «3i2/ - ^y 2 
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is positive on [7., 7+]. Since (3 > 0, it suffices to choose j± so that 
Q(h±) > 0. If 7+ > is sufficiently small, then Q( r y+) is positive since 
Q(0) — d\ + k\ x /2 > is, and if Q(y) is positive at y = —7 (which 
is condition ()3.14|) ). then it is positive for 7_ in a sufficiently small 
neighborhood of —7. □ 

5.2. Proof of Theorem 13.31 part II: verification of conditions 

(iv)— (v) of the general scheme. When x i n the definition of <f) e is 
fixed, the depends on e only, and to stress the dependence on e, in 
this subsection, we write L e instead of L^. In the proofs of the existence 
and the uniqueness theorem for the solution of problem (j4.1Uj) - (j4.11j) 
and estimates (j4.15j) and (j4.17j) . we use an additional property of r°, 
which was omitted from the general scheme in Subsection 14.11 as too 
technical (and by no means necessary), and which holds for and 
rjv,^- Fix we (0, 1/2). We denote by r e any function (in fact, a family 
of functions, parametrized by e G (0, 1)), which admits bound ()4.4|) 
with the constants C\ > 0, C 2 and M = M e _1 , where Co,Ci,M are 
independent of e G (0,1); in addition, r e must satisfy the following 
Lipschitz condition with weight: 

there exists a constant Mi > independent of e G (0, 1) and such that 
for any x, y G R + x R satisfying \xj — Vj\ < \ + e~ 1 ) w , j = 1, 2, 

\r e (x) - r e {y)\ < M x (|x| + e^e 1 ^. (5.8) 

Lemma 5.3. Functions r^^ and satisfy H4-4]) and uni- 
formly in N and e G (0, 1). 

Proof. Each of the functions and r^,^ is the sum of a linear function 
r or a piece- wise linear function rjv 5 for which ()5.8|) is evident since 
lu G (0, 1/2), and the function of the form 

f^x) = b 1>e (x 2 )x 1 + b 2 ,e(x 2 ) + b ,e, 

where b 0tt is bounded (hence, satisfies (j5.8J) ). and bj tt , j = 1,2, are 
constant outside the set x 2 G [— 1/e, 0], and have derivatives of order 
e on this set. By applying the Lagrange theorem, we obtain ()5.8J) for 
bi^ e (x 2 )xi + b 2 e (x 2 ). Thus, ()5.8|) has been proved. 

The proof of bound (|4.4j) in Subsection 15.11 is applicable to any r e , 
and the constants in (J4.4|) can be chosen the same for all e G (0, 1). □ 

The — L e is an elliptic operator in the half-space x\ > 0, which 
degenerates at the boundary x\ = 0, of the form 

-L e = -Mi - \x x dl + fi e (x)d 2 - ^P^dl (5.9) 
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where 9i,a,(3 are positive, and p e is an affine function of x\ 

/i e (x) = p efi {x 2 ) + x 1 n e>1 (x 2 ) > (5.10) 

whose coefficients are uniformly bounded, constant outside the segment 
[— 1/e, 0], and satisfy the global Lipschitz condition with parameter: 

\He,j(y) ~ Ve,j{z)\ < Ce\y - z\, (5.11) 
where C is independent of e G (0, 1). 

Theorem 5.4. Let —L e be an operator of the form \5. ty) . where 9\,a, j3 
are positive, and the real-valued function p t satisfies \5.1U\) and H5.ll)) . 
Let r € satisfy H4-4)) an d \5.8\) . with constants C, C , C%, M independent 
o/eG (0,1). 

Then there exist e > such that for all e G (0, eo) 

a) for any T > and any continuous g°, which exponentially decays at 
infinity, the problem 

(d T -L e + r e )f{x,r) = 0, rG(0,T), (5.12) 
f{x,r) = g°, (5.13) 

has the unique continuous solution, which exponentially decays at 
infinity; 

b) set Ao = 2M e _1 , o = n/2 + e; then there exists C\ and p > such 
that if \p\ < po then for all e G (0, eo) and X G £<t,A - 

1 1 exp(p(.)")(A + r°(-))(A - L t + r,)' 1 exp(-p(.)")ll < C v (5.14) 

Proof, a) Instead of the problem on a strip r G [0, T], we can con- 
sider the problem on r > 0, and look for the solution in the class of 
continuous functions which admit a bound 

\f(x, r)\ < Cexp[~W + r(X - (5.15) 

for some 5± > 0. By using the Laplace transform, we can reduce the 
existence and uniqueness theorem to the uniform invertibility of the 
family A — L t + r t , on the line 9ftA = Ao; the uniform invertibility follows 
from part b). Notice that in order to prove that the solution satisfies 
the boundary condition (J5.13j) in the strong sense, the invertibility for 
A G So- Ao and estimate (|4.15|) are needed. 

b) will be proved in the appendix. Notice that both (|4.15|) and (|4.17j) 
are implied by ()5.14|) . □ 
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5.3. Proof of Theorem 13. 71 For j = 2, . . . , n, fix 7+j > and 

7-,j < (5-16) 
Next, take a non- decreasing function Xj £ C°°(R) such that 

Xj(v) = 7-j, 2/<-l; (5-17) 



(5.18) 



and 



7-,i < Xj(y) < 7+j, VyGR. (5.19) 
Then, for any e G (0, 1), construct functions XjAv) = Xj( e u) an d 

and finally, set 



e {Xj). 



(j) e {x) = K 1X Xi + ^?Pj. 

Clearly, it satisfies ()4.3|) . Set 

5 := maxj/tn, max7 + max{— 7,- — 7_,}}; 

j>2 ' j>2 

it is positive since kxi > 0,7+j > and (|5.1(jjl holds. 

Lemma 5.5. For any e G (0, 1), there exists C e such that for all X\ > 
and x 2 G R, 

exp(-0 e (x''))|/oM;x,t)| < C e exp(-5(x 1 + (5.20) 
In particular, holds. 

The proof of this lemma is an evident modification of the proof of 
Lemma 15.51 in the case n = 2, and the rest of the proof of Theorem 
13.31 is a straightforward modification of constructions and arguments in 
Subsection 15.11 and Subsection 15.21 



5.4. Proof of Theorem l3~TTl Fix 7 < -7 and 7+ > 0, k,l 2 G R, 
and for e > 0, construct Xe an d ip e as in Sectional After that, define 



Set 



9 1 



0i 
O2 



X) = liXi + l 2 X 2 + Ipe(x-i). 



di 

2 



x 



Xi 

x 2 
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We have 

exp(-(f) e (x))Lexp((f) e (x)) = (6 1 - k 1 V, d 1 + I 1 } 

-(K 31 Xi + K 32 X 2 + K 33 X 3 )(d 3 + Xe(x 3 )) 

+^x 1 (d 1 + l 1 ) + —x 2 (d 2 + l 2 ) 

therefore 

= + ha 3 + (3 3lXl + (3 32 x 2 )d! 



2 J J 2 

J"=l,2 



+ ((/3n/i - Kn)xi - Ki 2 x 2 )9i + ((/3 22 / 2 - /t 22 )a: 2 - /t 2 ia;i)9 2 

+ [-«31^1 - ^32^2 - ^33^3 + («3 + Att^l + A^^Xe]^, 

and by denoting the columns of the matrix k, 11 as j = 1, 2, 
rw(*) = d,-{e\l l )-^(e^(x 3 ) + Xe (x 3 f) 

+[dx + (k 1 , - + K 31 Xe(x 3 ) - ^(ex' e (x 3 ) + Xe(x 3 ) 2 )]x 1 

+ [d 2 + (K 2 , I) - + K 32 Xe(x 3 ) ~ ^-(eXeixs) + Xe(x 3 ) 2 )}x2 

+ [d 3 + K, 33 Xe(x 3 )]x 3 . 

Now it is clear what the optimal choice of l\ and l 2 is. Indeed, it is 
necessary that in the formula for L^, the coefficients at dj, j = 1,2, 
must be non-negative, hence lj > Kjjjfijj- On the other hand, in the 
formula for r^^, it is better to have the coefficients at Xj, j = 1, 2, as 
large as possible. Equivalently, 

Kuh - + K21I2 and k 22 / 2 - ^22l2l + Kvi \ x 

should be as large as possible. But for j = 1,2, Kjjlj — f3jjl 2 /2 attains 
its maximum at lj = Kjj/j3jj, and since k 2 i is negative, k 2 \1 2 is attains 
its maximum on [k 22 /(3 22 , +oo) at l 2 = K 22 /j3 22 . Hence, the optimal 
choice is lj = Kjj//3jj,j = 1,2, and with this choice, the coefficients at 
Xj,j = 1, 2, become 

Kn K 2 \K 22 

< \ 031. 



and 



dl + OR~ + ~7i + ^lXefas) - -7T^X'e( X 3) + Xe(x 3 ) 2 ) 

tPll P22 * 



K 22 K\ 2 K\i 

d 2 + -5- + —5 + K 32 Xe{x 3 ) - —{eXeiXs) + Xe{X 3 ) ). 

^P22 Pll * 
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If e > and 7+ > are sufficiently small, and 7_ < —7 is sufficiently 
close to —7, then under conditions (|3.64|) - (j3.65|) . both coefficients are 
positive and bounded away from uniformly in x 3 . This property 
allows us to repeat all the proofs in Section |S] with small and evident 
changes. 

There are two subtle point is in the construction of local represen- 
tatives and local almost inverses in the appendix: when the set £7^ ■ 
intersects with the plane X\ = 0, we take a; ej in this plane (similarly 
in the case of the intersection with the plane Xi = 0), and if with the 
line Xi = X2 = 0, then on this line. In the first case, we freeze all 
the coefficients except for the ones in the expressions X\di and Xidf 
(similarly for the local representatives at points in the plane x?, — 0), 
and in the second case, all the coefficients except for the ones in the 
expressions x\d\ and xidf, I = 1,2. In the first case, the model family 
are operators on the half-line, and in the second case, on the quadrant 
{x\ > 0,X2 > 0}. The remaining details are straightforward modifica- 
tions of the corresponding steps in the proof for the family A\(2). 



Appendix A. Proof of Theorem 15.41 b) 

A.l. Reduction to construction of almost inverses. We start 
with the case p = 0, and in the end of the proof, indicate which changes 
need to be made for the case of small \p\. 

For each A G C satisfying |A| > e -1 , introduce the family of Banach 
(in fact, Hilbert) spaces H\ C H := L 2 (R+ X R), which consists of 
functions with the finite norm 1 1 • \ \ H i defined by 

2 

:= IKkl + IADwII' + ll^lP + llxi^lP + ^dl^lP + llxi^H 2 ), 

j=0 

where || • || is the norm in H. Notice that the embedding H\ C H is 
continuous, and H\ is independent of A as a topological space; but it 
is convenient to derive estimates for the auxiliary operators below in 
terms of the norm depending on the parameter. 

To prove the invertibility of — L e +r e +X and estimate ()4.15|) . it suffices 
to construct left and right almost inverses, R l Xe and R r x t - These are 
operators which are uniformly bounded 

ii^yi^ < c- (A.i) 

\\R r xJ\ H ^Hl < C, (A.2) 
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where C is independent of e G (0, 1) and A G So-,a , an d satisfy the 
approximate equalities in the definition of the inverse: 

R[ e (-L e + r e + X) = I + T{ >£ , (A.3) 

(-L e + r e + X)Rl !e = I + T r Xt , (A.4) 

where 

\\TU\h^h -»• 0, ase^O, (A.5) 

IKMh^hI - 0, ase^O, (A.6) 

uniformly w.r.t. A G Xo-,a - Notice that (|A.4|) and (|A.5|) imply that 
_R Ae maps i7 into the domain of — L e + r e . Suppose that the almost 
inverses are constructed, and e is so small that the norms of operators 
T{ e and T{ e are less than 1/2. Then the RHS in (TOJ) and (TOJ) are 
invertible bounded operators, and hence — L t + r e + A has the bounded 
inverse. Moreover, the inverse R\ e (I + Tx^-r)' 1 is a bounded operator 
from H to H\, uniformly in A, which implies (|5.14|) . 

Thus, it remains to construct left and right almost inverses. The 
standard technique consists of the construction of 

(i) an appropriate partition of unity, 

(ii) local representative of the operator (localization of the operator), 
that is, the freezing of coefficients w.r.t. to all the variables or 
some of them at some points, 

(iii) the inverses to the local representatives, and finally, 

(iv) global almost inverses by using the partition of unity and the local 
inverses. 

Notice that Ti a> \ depends on e, and |A| — > +oo as e — > 0, uniformly in 
A G £<t,Ao ! this property will be used systematically in the constructions 
below. 

A. 2. Partition of unity. We use the universal construction due to 
Hormander (1985), which is based on the definition of a slowly varying 
metric; for the variant of the construction for operators with parameter, 
see Levendorskii (1993). Fix u G (0, 1/2), and for each e G (0, 1), define 
a function (y) e = (|e|~ 2 + y 2 ) 1 ^ 2 on R, and the Riemann metric G e on 
on R 2 by 

G eVB (z) = (m>r 2 i2ii 2 + (x 2 ): 2u \z 2 \ 2 . 

It is straightforward to check that the derivatives of the function (y)" 
w.r.t. y are bounded uniformly in y G R 2 and e G (0, 1), therefore (see 
Levendorskii (1993)) there exist c, ci, C\ > such that if G e - X (x— z) < c, 
then 

ciG e]Z (w) < G e . x (w) < C 1 G e , z (w), V w G R 2 . (A.7) 



34 



S. LEVENDORSKII 



Condition (|A.7|) means that the metric G e is slowly varying (uniformly 
w.r.t. e G (0, 1)), and therefore there exist positive constants C 2 , C 3 , C 4 , 
points x 6 '- 7 G R 2 , and non-negative functions v e j G C°°(R 2 ),j = 
1,2,..., such that 

(i) sets U e j = supp u e> j, j — 1,2, ... , cover R 2 ; 

(ii) the multiplicity of the covering {U e j}j>i is bounded uniformly 
w.r.t. e G (0, 1), and j — 1,2, ... , that is, for each j, the number 
of k for which U e j and £7^ intersect is bounded by C 2 ; 

(iii) x £J ' G I7 e j; 

(iv) the diameter of U e j is not greater than C*2 ( |e| — 1 + |x e '-'|) w ; 

(v) for all multi-indices s, e G (0, 1), and j > 1, 

\d s u £>j (x)\<C s (x^):^, Vx, (A.8) 

where \s\ = si + s 2 , and the constants C s are independent of 
e G (0, 1), and j > 1; 

(vi) for all e G (0, 1), and all x G R 2 , 

5> eJ (x) = l. (A.9) 

For details of the construction, see Hormander (1985) and Levendorskii 
(1993). In the construction, one can choose the points and functions 
so that 

(vii) either x eJ is on the line x\ = 0, or x eJ and U t j are outside the 
strip |sei| < C2 1 e | where c 2 > is independent of e G (0, 1) and j. 

From now on, we assume that (i)-(vii) hold, and we consider only the 
points in the half-plane x\ > 0; now the v e j are functions defined in the 
same half-plane, and (jA.9|) holds for x from this half-plane. We divide 
the set of the points x e ^ in the half-plane x\ > into two subsets: 
j G Jo, if is on the boundary X\ = 0, and j G J+, if x e 'i belongs to 
the open half-plane X\ > 0. 
Construct 

and set U}^ = supp U t j. Similarly, starting with v\ - and U}p construct 
vi - and U^j Then v\ - and U l e j, I = 1, 2, satisfy the same conditions as 
u e> j and U t j (with different constants) but (|A.9[) . In addition, 

vljix) = 1, VxeU £tj , (A.10) 
vl 3 {x) = 1, Vie^. (A.ll) 
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A. 3. Local inverses to non-degenerate local representatives. 

For j £ J+, denote by A Xjt j the operator — L t + r e + A with the coeffi- 
cients freezed at x €,3 \ 

Ax,e,j = -Old, - -x?dl + fi e (x^)d 2 - d 2 2 + r e (x^) + A. 

The function 

1 cm 4- Br 6 '-* 

o w (0 = -Witi + 2 X ?£ + + 2 1 £2 + r e (s eJ ') + A 

is the symbol of the operator A Xjt ,j, that is, 

Ax, e ,j = a\,e,j{D) := ^ r ~ 1 a Aie)i (0^ r , 

where JF is the Fourier transform. Equivalently, for a sufficiently regular 
function u, 

a x , e>j (D)u(x) = (2n)- 2 [ e«a w (£)«(£)d£, 

where 

is the Fourier transform of it. 

Lemma A.l. There exist c > and eo > stzc/i that for all e £ (0, eo), 
A £ S^Aq, j £ J+ and £ £ R 2 ; 

|aw(0l>c(|A| + |x«1+x^| 2 ). (A.12) 

Proof. Due to our choice of cr, 

3feA > -ei|9A|, VA£S CTiAo , (A.13) 

where ej — > as e — > 0. Hence, if C > is fixed, then for all A £ S^,,, 
j £ J + and £ £ R 2 , satisfying 

|A|<C(|^'|+^|^| 2 ), (A.14) 

we have 

1 , . rv 4- /9r eJ 

and for sufficiently small e > 0, (|A.12|) obtains. On the other hand, if 

\\\>C(\x^\+x\' j \tf), (A.15) 

then \ax, e ,j(0 ~ -^ll^l -1 — ► 1, as (7 — > +00, and hence if C is large, and 
A £ S CT;Ao , j £ J + and £ £ R 2 satisfy condition (jA.15j) . then ()A.12|) 
holds as well. □ 
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In view of (1A.12J) . we can define pseudo-differential operators (PDO) 
R\, e ,j,j e J+, by 

Rx,ej = a x ^{D)- 1 := T^ax^ii)- 1 ^. 

By using the Fourier transform, we immediately obtain that Rx,e,j is 
the inverse to A\ ie> f 

Ax,ejR\,e,j = I, R\e,j^\,e,j = I- (A.16) 

A. 4. Local inverses to degenerate local representatives. For j G 

Jo, define 

A\, e ,j = -Bid, - \x x d\ - ^±I^± d l + n {x^) + A, 

and introduce the Hilbert space of functions on the half-plane 
x\ > 0, with the finite norm || • 1 1 A,e,j defined by 



?! 



ll e j = llft«l|V||x a ^tt|| a +^(||^u||V||xi^u|| a )+||(|rr^|+|A|H| a l 

s=0 



where || • || is the norm in H = L 2 (R + X R). 

Lemma A. 2. There exist eo > and C such that for each A G £ CT ,Ao; 
the operator Ax, e ,j : Hx, e ,j ^ H is invertible, and the norm of the 
inverse, call it R\, e ,j, is bounded by C . 

Proof. This a very special case of general results for degenerate elliptic 
operators in Levendorskii (1993). We recall the scheme of the proof. 
First, we make the (partial) Fourier transform T w.r.t. to X2, and 
obtain 

A x>eJ = f- l A x ^ 2 )F, 

where 

x „ d Id 2 a + By ~ , , ... 

is the family of operators on R + parametrized by £ 2 £ R an d A G S CT Ao . 
Introduce the Hilbert space Hx, e ,j,& °f functions on R + , with the finite 
norm || • 1 1 A,e,j,^ 2 defined by 

IMIL,, 6 = IKI! 2 + li^'ii 2 + ll(i + v)?«\\ 2 + IKbSfl + WHi 2 , 

where || • || is the norm in L 2 (R + ). By using an appropriate partition 
on unity on R + , which depends on (A, £2), one shows that if eo is 
sufficiently small then for all e G (0, eo), A G £<t,a j an d £2 G R, operator 

Aa,6j(6) : #A,e,j,6 -> ^2(R+) 
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is invertible, and the norm of the inverse, Rx^jfa), is bounded by a 
constant which is independent of A, e, £ 2 - 
It remains to define 

□ 

A. 5. Construction of almost inverses. Set 

R \e = ^2 U e,3 R \e,3 U lj> 
J>1 

R \t = ^2 U eJ R \eJ u eJ- 

We check that if eo > is sufficiently small then the conditions (|A.2|) . 
flA.4|) and (|A.5|) hold for the family of almost right inverses R r Xe ] the 
verification of their analogs for almost left inverses R l x e is similar. 
Notice that 

1) the multiplicity of the covering {U^j}j>i is finite; 

2) the norms are defined in terms of integration of a function and its 
derivatives multiplied by weight functions; 

3) the derivatives of functions V s • satisfy estimates ()A.8|) . and 

4) (|A~9jl . (IA~10|) and ljXlT|) hold'. 

Therefore, it suffices to prove that 

a) there exist eo > and C such that for all e G (0, eo), A G £„- \ , and 
j > 1, 

lk^V,,Hif^<^ (A.17) 

b) let eo — * 0; then uniformly in e G (0, eo), A G S CT ,Ao) an d j > 1, 

| |^-(-L 6 + r e + A - .1 A ,.,!<) I - 0; , (A.18) 

and 

- A Xi€ y ed )\\ H i_> H - 0. (A.19) 

Indeed, (|A.2|) follows from (|A.17|) since the multiplicity of the covering 
{U}j}j>i, call it M, is finite: for any u G H, 

\\Rx >e , r u\\ H i < MsupWvljRx^jVejuWH 
j 

< MC sup 1 1 U e jU 1 1 H 
j 

= MC\\u\\ H . 
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Further, by using IjA.llj) . we can write 

(-L e + r £ + X)R x ,e,r = ^2^. r \,\.,.jf>.\.,j^j (A.20) 

i>i 

i>i 

By construction, A\ )€) jR\ )€) j = I, therefore due to (jA.9|) and (jA.lOj) . 
the first sum on the RHS of (|A.20|) is identity: 

j>i i>i 

Since the multiplicity of the covering {U^Aj>\ is finite, we conclude 
from (lA~P7l) . dUD) and (lA~l9|) . that the norms of the last two sums 
on the RHS of (jA.20|) . as operators in H, tends to as eo — > 0, uniformly 
in e G (0,eo),A G Xo-,a - This proves (jA.4|) and ()A.5|) . and finishes the 
proof of part b) of Theorem 15.41 

A.6. Proof of estimates ()A.17jh ()A.18j) and ()A.19j) . It suffices to 
use the following simple observations: 

1. Due to the choice uj G (0, 1/2), for j G Jo, f|5.8|) gives 

||(r e (x)-r £ (^))^.||^ if ^0 

as e — > 0, uniformly in j; 

2. Since function /z,, is an affine function of x\ with coefficients 
uniformly bounded w.r.t. x 2 , e G (0, 1) and the other parameters, 
and the coefficients satisfy (j5.11|) . we conclude that 

as e — > 0, uniformly in j; 

3. For j G J_|_, an error of freezing of any coefficient at x 6 ^ , on the 
set U^j produces a relative error of order e"; 

4. Any derivative of any v e j and i/' • is of order e^, and hence each 
time we calculate such a derivative, the term, in which the derivative 
enters, has a small norm. 

A. 7. Proof of (J4.17|) for small p. Set L e ,p = ^v{p{ x 2))L e exp(— p{x 2 }) 
Direct calculations show that L e>p — > L e as p — > 0, in the operator norm, 
uniformly in e, therefore for small p, the almost inverses to A — L e + r e 
can be used to construct the inverses to A — L e ,p + r e , and the rest of 
the proof of ()4.17|) remains the same. 
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